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Abstract 



In the previous work (J. Geom. Phys. 39 (2001) 50-61), the closed loop sohtons in a plane, 
i.e., loops whose curvatures obey the modified Korteweg-de Vries equations, were investigated for 
the case related to algebraic curves with genera one and two. This article is a generalization of 
the previous article to those of hyperelliptic curves with general genera. It was proved that the 
tangential angle of loop soliton is expressed by the Weierstrass hyperelliptic al function for a given 
hyperelliptic curve = f{x) with genus g. 



This article is on loop solitons related to hyperelliptic curves with higher genera as an extension 
of the previous report [Ma5]. 

In [Ma2], I proposed a problem of a quantized elastica (ideal thin elastic curve), or statistical 
mechanics of elasticas, which is a model of a large polymer in a plane, such as DNA at finite 
temperature. When a position of the elastica in a complex plane C is denoted hy Z : ^ C, the 
partition function of elastica is given by 



where /? is the inverse of temperature, DZ is a certain functional measure and E[Z] is the Euler- 
Bernoulli functional energy. 



for the arclength s, i.e., induced metric of the curve, and its curvature k. 

As shown in [Ma2-4] , the partition function is completely determined by the orbits of the modified 
Korteweg-de Vries (MKdV) hierarchical equations. As a curve obeying the MKdV equation is 
known as a loop soliton due to [KIW, I], the quantized elastica problem is a realization of the loop 
solitons. It is worth while noting that even though the soliton theory is studied in field of physics. 
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there are not so many examples that sohton equation is connected with a physical model including 
its multi-soliton solutions. 

In this article, we will consider hyperelliptic solutions of the loop solitons or excited states of a 
quantized elastica. Theorem 3-2 is our main theorem of this article. There we give explicit solutions 
of closed loop solitons in a plane related to hyperelliptic curves with general genera. 

We will base on the result of [Ma7]; there we show hyperelliptic solutions of the MKdV equation 
in terms of the theories of the hyperelliptic functions which were developed in nineteenth century 
[Bal-3, Kl] and are recently re-evaluated [BEL1,2, Ma5, 6]. Following the idea mentioned in 
the discussion in [Ma5], we extend the results in [Ma5] to the case of general genus in terms of 
Weierstrass's hyperelliptic al functions [Ba2 p. 34, W]. 

As the elastica problem has a deep history [L, Tl, 2], I believe that one of its reasons is its 
naturalness. In the derivations of the solutions, it turns out that the elastica problem is very natural 
even from a mathematical viewpoint. In Remark 3-3, we will give comments on its naturalness. 

Further as I mentioned in [Mai], the elastica problem is closely related to automorphic function 
theory even though the solutions are constructed in Abelian variety. In fact the Euler-Bernoulli 
energy functional can be expressed by the Schwarz derivative, 

E[Z]= ^ ds{Z,s}sB 

In §4, we comment on its relation to automorphic functions. 

§2. Differentials of a Hyperelliptic Curve 

In this section, we will review the hyperelliptic functions following [O, BEL2, Bal-3] without 
explanations and proofs. 

We denote the set of complex number by C and the set of integers by Z. 

Convention 2-1. We deal with a hyperelliptic curve Xg of genus g {g > 0) given by the affine 
equation, 

y' = fix) 

= A2g+lX^»+^ + \2gX^3 + ■■■ + \2X^ + \lX + Ao 

= P{x)Q{x), 

where 

f(x) = {x~ bi){x - 62) • • • (a: - b2g){x - 629+1), 
Q{x) = (x - ci)(x - C2) • • • (x - Cg){x - c), 
P{x) = (x - ai){x - 02) • • • (x - Qg), 

'^2g+i = 1; o-nd \j 's, aj 's, bj 's, cj 's and c are complex values. 

Definition 2-2 [Bal,2, BEL1,2 6]. 

For a point {xi,yi) G Xg, we define the following quantities. 
(1) Let us denote the homology of a hyperelliptic curve Xg by, 

g g 

Hi(X„Z) = 0Za,e0Z^,-, 
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where these intersections are given as [ai.aj] = 0, [l3i,Pj] = and [ai,(3j] 
(2) The unnormalized differentials of the first kind are defined by, 



7 (i) XidXi 



(3) The unnormalized period matrices are defined by, 



du] 



(a) 



du) 



(a) 



U3 :- 



Ui" 



(4) The normalized period matrices are given by. 



' * [ du{'' ■ ■ ■ du. 



(»)] 



r := a; u , u :- 



^9 
T 



(5) The unnormalized differentials of the second kind are defined by, 

xj^'^dxi 



""l •— ) ""■2 7t : 



du^: 



2yi 



and dr(') := {dr^\dr^\- ■ ■ , dr^g^). 



where A is 2g x g matrix defined by 



/O A3 2A4 3A5 
As 2A6 



A 



A7 





{9 - gXg+2 {9 + l)A3+3 

(g - 2)Ag+2 {g - l)Xg+3 g>^g+4 
{g - 3)Ag+3 {g - 2)Ag+4 {g - l)Ag+5 



A2g-2 





2A 



A 



2g-l 



2s+l 



3A2g+l 





(25 - 3)A2g-i (2<7 - 2)A23 (25 - l)A2g+l \ 

(25 - 4)A2g {2g - 3)A2g+i 
{2g - 5)A2g+i 








(6) The complete hyperelliptic integral matrices of the second kind are defined by. 



ri :-- 



dr. 



(a) 



^3, 



V 



dr. 



(a) 



^3. 



UJ 



(J 



(7) By defining the Abel map for g-th symmetric product of the curve Xg, 

u:Sym^(X,)^C^ (^km)i=i,--- ,9) ■=Y.f '^'^^k ^ k = l,---,^, (2-11) 

the Jacohi variety Jg are defined as complex torus, 

Jg := CVA. (2-12) 

Here A is a lattice generated by uj. 



Definition 2-3. 

The coordinate in O for points {Qi = {xi,yi) \ i = 1, ■ ■ ■ , g } of the curve y^ = f{x) is given 



by, 



3 Axi,yi) 9 

u.:=J2 duf, duj = ^duf. (2-13) 

i=i i=i 

(1) Using the coordinate Uj, a functions, which is a holomorphic function over , is defined 
by [Ba2, p.336,p350, BELl, Kl], 



a{u) = a{u;Xg) := 7exp(— ^ * wrf'u' ^u)-d 
where ^ is a certain constant factor. 



6" 
6' 



{z;t) := ^ exp 



27rv^ \-Xn + a)T(n + a) + \n + a){z + h) 



for g-dimensional vectors a and b, and 



S' :-- 



2 2 



6" 



1 1 
2 



(2) Hyperelliptic p-function is defined by [Bal, Ba2, Kl], 

(Pij{u) := - logc7('u), 

OUiOUj 

and Hyperelliptic Q function is defined by, 

(3) Weierstrass hyperelliptic air-function is defined by [Ba2 p. 340, W] 



a\r{u) :=7V^(M> 

where 7' is a certain constant, 

F{x) := (a; — xi) ■ ■ ■ {x — Xg) 

=: 7gx3 -I- 7g_ix^"-^ H h 7o- 



(2-14) 



(2-15) 



(2-16) 



(2-17) 
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(2-20) 
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Proposition 2-4. 

(1) pgi {i = 1, - ■ ■ ,g) is elementary symmetric functions of {xi,X2, ■ ■ ■ ,Xg} [Bal, 2, 3, BELl], 



I.e., 



F{X) = X» - ^ 



i=l 



(2) Cj{u) is expressed by [BEL2 p.33-35] 



k=l 



dri - -detA 



(2-21) 



(2-22) 



where 



A 



ra+l 





ei 


-1 








\ 




262 


-ei 


1 









(n 


- 2)e„_2 


-en-3 


6n— 4 


• ±1 






(n 


- l)e„-i 


-e„-2 


en-3 ■ 


• ±ei 


Tl 




\{n 


ne„ 


-e„-i 


en-2 • 


• ±62 


T6l 


±1 


+ l)c^n+i 




dn-1 • 


• ±^3 




±diJ 



(2-23) 



and 6i := anc/ := pg,g,g+i_i := dpg^g+i-i/dug. 

As we show later, (2-22) is very important relation in our quantized elastica, which was found 
by Buchstaber, Enolskii and Leykin [BEL2] (Appendix); according to [BEL2] Baker in [Bal] gave 
a wrong relation. We note that F{x) is a generator of the elementary symmetric functions and the 
matrix An (2-23) contains the matrix of the Newton formula as its minor matrix. In fact in the 
derivation of (2-23), A^ plays the role which connects the elementary and power sum symmetric 
functions. 

Definition 2-5. 

(1) A polynomial associated with F{x) is introduced by 

F{x) 



TTi{x) := 



Xi,g-lX^ ^ + Xi,g-2X^ ^ H h Xi,l + Xi,0, 



(2-24) 



where Xi,g-i = 1; Xi,g-2 = {x^ + 
(2) We will introduce g x g-matrices 



W 



/Xi,o Xi,i 

X2,0 X2,l 



+ Xa 



Xl,g-1\ 
X2,g-1 



y 



^Xg,0 Xg,l ■■■ Xg,g-\' 
/F'(xi) 

V 



F'{X2) 



Xi, and so on. 

\ 

\ 

F'{Xg)j 



y2 



Vg 



(2-25) 



(3) 



where F'{x) := dF{x)/dx. 



( 1 

7s- 1 1 

7g-2 79-1 1 



M :-- 



\ 



7i 72 







7.9-1 1/ 



K 



g-1 g-2 



1\ 
1 



(2-26) 



(4) r/ie coordinate in is introduced by u^^^ := "P^u, where Vr is defined by its inverse matrix, 



5-1 



2 



1 56^ 

1 {g- l)h 



9 1^59-1 5s-i\ 

g-lj"r 

g-2 J 





yo 








1 




(2-27) 



1 / 



(5) For a polynomial g{X) = gnX'^ + ■ ■ ■ + go, we introduce the Dj operator 



(2-28) 



Lemma 2-6. 

(1) The inverse matrix of W is given by = T~^V , where V is Vandermonde matrix, 

/I 1 ... 1 X 



V 



Xi X2 

™2 ™2 
Xi X2 



X 



\xl 



9-1 



(2-29) 



(2) Let du, := d/dui, := d/dxi and d^^ := d/duf , 



(3) KM = W and 



3(0 



t-p-1 



Wij=Xi,j-i = [Dj{F{X))]x=.^. 



(2-30) 



(2-31) 



Proof. (1) is obvious by using the properties of the Vandermonde matrix. In (2), we must pay 
attention the fixed parameters for the partial differential. By comparing c?x* and the chain relation 

/ Fix) \ 

of , we obtain the matrix representation (2-30) [Ma7]. Prom the relation ( - — — ) • {x — Xi) = 



F{x), we have 

Xi,j — Jj+l + XiXi,j+l- 

Then we obtain the relations (2-31). ■ 

We note that the formulae (2-31) and (2-32) are very important to prove (2-23). 



(2-32) 



Proposition 2-7 [BEL2 p. 11 ]. 

The Legendre relation is given by 

^Lo'ri" - *a;"r?' = 2Tr V^Ig, (2-33) 

where Ig is the g x g-unit matrix. 



3. Loop Solitons 



In this section, we will deal with a real curve in a plane in the category of differential geometry. 
Let us consider a smooth immersion of a circle into the two dimensional Euclidean space 
~ C or -f- {oo} ~ CP^. The immersed real curve C is characterized by the affine coordinate 

{X^ (s) , (s)) around the origin. Here s is a parameter of and is, now, chosen as the arclength 

so that ds'^ = {dX^)^ + {dX'^y. We will also use the complex expression, 

Z{s) := X\s) + V^X^. (3-1) 
Then by letting dg := d/ds, \dsZ{s)\ = 1 and the curvature of C is given by 

k{s) := ^d,logd,Z{s). (3-2) 



As mentioned in the introduction, loop soliton is identified with a quantized elastica [Ma2]. Thus 
we will sometimes call it quantized elastica or simply elastica hereafter. 

Definition 3-1. 

(1) A one parameter family of curves {Ct} for real parameter t is called a loop soliton, if 
its curvature obeys the MKdV equation; for q := k/2, 

dtq + Qq^d.q + = 0, (3-3) 

where dt ■= d/dt. 

(2) The energy of elastica is given by 

E[Z] := jdsk^, (3-4) 
which can be expressed by the Schwarz derivative, 

E[Z]= ^ds{Z,s}sD, (3-5) 

where 



Here we will give our main theorem. 
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Theorem 3-2. 

Let the configuration of the x- components {xi, • • • , Xg) of the affine coordinates of the hyperelliptic 
curves Sym^(Xg) satisfy, 

\Fibr)\=ro, (3-7) 

where rp is a positive number. For such (xi, yi), • • • , {xg,yg), we have u := u((xi, yi), • • • , {xg,yg)) 
due to (2-11). 

(1) By setting s = Ug/ro and t = Ug-i + {\2g-1 + br)ug, 

:=F(6,), or = 1, 

completely characterizes the loop soliton. 

(2) The shape of elastica is given by, 




Remark 3-3. 

(1) If one prefers more proper expression for the branch point {br,0), he may use Uj~'' in (2-27) 
and then find similar results. As the expression is essentially the same as the above one due 
to (2-30), we will investigate the above one. 

(2) The condition (3-7) is essential. Due to the condition, any configurations of u, the tangential 
angle ^ = \og{dug Z) / does not contain imaginary part. Hence the arclength locally 
does not change [Mal-4]. As Goldstein and Petrich showed that the isometric deformation 
of space curve in a plane gives the MKdV equation [GP], this condition and = 
recover the MKdV equation in general. 

(3) The tangential vector, du^Z = F{br) oc al^, consists only of Xj's, which can be regarded 
as a two- fold coordinate of Sym^(CP"^). Each Xi G CP^ appears when we construct the 
hyperelliptic curve Xg using two CP^ with g -\- 1 cuts. 

(4) Due to the configuration of (3-7), there is a trivial action of the U(l)-group, which exhibits 
the translation symmetry of the elastica. When we begin with this symmetry and isometric 
deformation, we reproduces the MKdV hierarchy [Ma2,3]. 

(5) The condition (3-7) should be regarded as a reality condition. Elastica problem is a real 
analytic problem. In the primitive sense, the complex analysis is more complex than the 
real analysis but from deeper viewpoint, their standpoints are reversed. In fact, due to 
the condition, we must investigate all possible contours in the complex curve Xg. In other 
words, from the point of view of real analysis, as long as we have insufficient knowledge of 
the condition, it is not the end of the study of the quantized elastica problem. I suppose 
that this difficulty is similar to that of real analytic Eisenstein series [S]. 

(6) The condition (3-7) is satisfied if all Xi are in a circle centrizing at 6^. Then symmetric 
configuration of xi, • • • , Xg determines a point of a shape of the loop soliton. In other 
words, the dynamics of the elastica is translated to symmetric system of (/-particles in S^. 
Dynamics of symmetric particles in a circle might be familiar with researches of quantum 
integrable system [KBI]. 

When we consider the discrete configurations of x's, they give the discrete time develop- 
ment of the piecewise linear curves. This must be related to the discrete integrable system. 

From the definition. Theorem 3-2 can be proved by the next proposition, which was shown in 
[Ma7] . We will give a sketch of the proof of [Ma7] , whose techniques essentially appeared in [Ba3] . 
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(3-8) 



Proposition 3-4. 

By letting 



jjL^^^ obeys the modified KdV equation, 



(3-10) 



(3-11) 



Proof. This is proved by [Ma7] . We will give a sketch of the proof. Prom the definition, we have 



_d_ 

dUr 



logF(6.) = ^- 



Let dX° is boundary when Xg is embedded in a upper half plane H. By estimation of 



exo {x - br)F{xY 



dx = 0, 



and counting its residues, we obtain 



d 



fix) 



9 1 

^ F'{xk) [dx \ {x-br)F^{x 



= Mg + br + "^Pgg- 



J X=Xk 



Further we have 

Vk 



E 



k 



{x - Xk)F'{xk) 



= E 



'^VkVi 



kik^ii^ ~ Xk)ixk - xi)F'{xk)F'{xi) ' ^ (x - Xk)^F'(xky 



Using them we have the relation. I 
We note that the formal power sires 



(r) — 



I d_ 



logF{br) 



1 



1 ^ F' 



00 g 



{Xi){Xi - br) 



i=l ^ '^v ' ' ' V 

is resemble to the generator of the power sum symmetric functions. 



Vi 5. 

{Xi)br bi- 



(3-12) 



(3-13) 



(3-14) 



(3-15) 



(3-16) 



As we proved Proposition 3-4, we give two corollaries, which are shown by direct computations. 
Corollary 3-5 gives local properties of the elastica and Corollary 3-6 is associated with its global 
properties. 

Corollary 3-5. 

(1) The shape of elastica is given by 



Z(-) = - +Y.^r dn-i + 2 H det 



(3-17) 



(2) The Schwarz derivative of Z with respect to Ug 

{ZW,nJsD = 4p9g + 2A2.g + 2hr. (3-18) 

(3) The root square of the tangential vector d^^Z^'^^ = aV is a solution of the Dirac equation 
or Frenet-Serret equation [Mai], 

= 0. (3-19) 






Here we comment on (3-18) and (3-19). First we note that the solution of the Dirac equation 
consists of al-functions due to (2-19) and (3-8). Second it is noted that from the definitions (3-6) 
and (3-10), (3-18) agrees with the Miura transformation, 

^ir)2 ^ ^Q^^fj^ir) ^ 2pgg + X^g + K, (3-20) 

because the left hand side consists of the solutions of the MKdV equation (3-11) whereas the right 
hand side obeys the KdV equation [BELl, Ma6]. Further it is obvious that (3-19) has the same 
data as the Miura transformation (3-20) by operating the Dirac operator twice [Ma7]. On the other 
hand (3-19) can be expressed by 

Jcl. \ / n ..(rW ( Jd,. \ 

(3-21) 

Here we can recognize that the left hand side is an operation in analytic category while right hand 
side is an operation as an endmorphims in a commutative algebra. This relation is essential in 
study of P-module, due to the statements in p. 12-13 in [Bj]. 

Corollary 3-6. 

(1) The winding number of elastica can be computed for a given path by the integration, 

w := ^ y du^(p{ug)dug. (3-22) 

(2) The closed condition of elastica, 

jd^^Z^'-^dug = 0, (3-23) 

consists of the conditions 

6(^)% + f^6WVi=0, (3-24) 

i=l 

using the hyperelliptic integral (2-5) and (2-6) 

(3) 

{Z^'-\ug}sudug = -AdCg + 2{X2g + hr)dug. (3-25) 

(4) 

/ {Z^^\ug}si,dug = -An'^g + 2{X2g + brW^g. (3-26) 

i {Z^'-\ug}sT,dug = -4.n',g + 2{X2g + brW^g. (3-27) 



4. Discussion 



As mentioned in [Mai], our problem is resemble to Poincare, Klein and Schwarz theories of the 
automorphic function over the half plane or the Poincare disk. 

First the resemblance is obtained through the conformal field theory. Since is homotopically 
equivalent with C — {0}, most results of the conformal field theory are obtained by investigation 
of dynamics of functions over in a Riemann sphere CP^ [Ka]. The situation in (3-7) is very 
resemble to that of the conformal field theory but is in the higher genus Riemannian surface. 
However we can also perform the Fourier transformation or localization around the ramified point 
{br,0). Then it should be noted that on the injective maps, 

there are nontrivial actions of the fundamental group (3-23). By letting 

K := j {2pgg + \2g + K) e^-V^-./'^o (4.2) 

we have the relation of Virasoro algebra [Be, Mai], 

[L„, Lm] = {n- m)Ln+m + So,n+mn{n^ - 1). (4-3) 

Here we remark that in the conformal field theory, the surface obtained by means of the Vertex 
operator acting a Riemannian sphere is called a Riemannian surface with genus g but as in showed 
in [Ma6], such a Riemannian surface, at least of the case of the hyperelliptic curve, is very far 
from our Riemannian surfaces. The Riemannian surface obtained by the Vertex operator from 
the Riemannian sphere is, in fact, topologically genus g surface but is very special (semistable) 
degenerate curve; the theory over such a curve should be regarded as a theory on Riemannian 
sphere [AE]. Recently some of physicists might regard that topological aspect in field theory is 
the most important and they deal only with degenerate curves. However, at least, in low energy 
physics related to our lives, we need finer topology. In fact, the shape of classical clastica, which 
was studied by Euler as a classical field theory and is determined by curvature, leads us to very 
fruitful physics and mathematics. Further in general, physical phenomenon are not in a framework 
of complex analysis. Even though some objects in the category of the complex analysis is classified 
by topological objects, it is not all in physics. For example, Euler equation for complete fluid 
dynamics in three-dimensional space cannot be expressed by complex analysis. Quantized clastica 
problem should be also considered in the real analytic category as Euler did. Of course, topological 
aspect is still important but, I believe, is not a goal for quantitative science. 

Next I will comment on interesting relations, which also looks connected with the automorphic 
functions. Z*^*") is roughly equal to (^-functions due to (3-17) whereas the incomplete "energy 
integral" , 

p {Z^^\Ug}si,dUg, (4-4) 

is also expressed by the C, functions from (3-18). It implies that there is a similarity between 
the configuration and energy of elasticas. Further the matrix (2-23) is essentially the same as 
the Newton formula which connects the elementary and power sum symmetric functions. Thus 
^, (the configurations and energy from above view points), is essentially expressed by the power 
sum symmetric function, while the main part of fi, or a half curvature of elastica, is the same as a 
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generator of the power sum symmetric function due to (3-16). It is expected that there might be 
hidden symmetries. These facts might remind us of rcplicability of automorphic function if they 
were related to automorphic functions [Mc]. In fact, for the case of elhptic function (g = 1) case 
helping with the notations in [Ma5], ((5„p(m))^ = 4p^ — g2p + gs = 4(p — ei)(p — es){p — 64), 
p{uji) = e, ), we have interesting formulae, 

d^Z^'^^u + cJa) = ^{Z(«)(u),u}sD - ^ei, (4-5) 



/u 1 -, pu pu 

d^__exp(--y^ [{Z^-\u"),u"}sB 

- {u" - LOa), «"}sd] du"du'^ , 



(4-6) 



for a = 1, 2, 3. In the integrations, we should note the effects from the initial points. These 
equations (4-5) and (4-6) can be extended to higher genus. For the case of (4-6) we can do by using 
the relation al^ functions and a function [Ba3]. 

Further from the point of view of theory of the symmetric function, one might wish to regard 
Xi as an eigenvalue of some matrix, X = diag(a;i, 0:2, ■ • • ,Xg). The generator of the elementary 
symmetric function is expressed by, 

F{x) = det{X -xl), (4-7) 

and that of the power symmetric function is expressed by 

00 

G(x) = ^tr(A^")x'^. (4-8) 

n=l 

They are closely related to (2-20) and (3-16). As comparison with physics, some purposes in 
mathematics are classifications and determination of the relations among classified objects. The 
classification should be characterized by discrete quantities and these discrete quantities should 
sometimes preserve when we take a certain limit. Thus it might be natural to consider degenerate 
curves in a certain sense. The degenerate curve = P{x)'^x, which was dealt with in [Ma6] 
and is associated with the soliton solutions and algebra of vertex operators, is also expressed by 

— a;(det(a;J — P))"^ = by letting V := diag(ai,a2, • • • ,ag). These matrices V and X for the 
degenerate curve are the same rank. Hence it is natural that one consider a graded algebra A 
{Ag C Ag+i, A = ^gAg) geucratcd by generators {X, V and so on) such that there is a ^-module 
Mg satisfying XMg = XMg and VMg = VMg. If we regard Mg as a representation of shape of 
elastica, we might able to deal with family of elastica related to hyperelliptic curves with different 
genera. The Schwarz derivative, which plays important roles in theory of automorphic functions 
[Mc] and is invariant for PSL(2, C), naturally appears in our elastica problem [Ma4]. Even though 
PSL(2,Z) is far from our situation in this stage, I hope that our study might reveal some relations 
between elastica problems and automorphic function theory [Mc], if exists [Mai]. 

Finally we mention our future study. The quantized elastica in a plane was extended to that in 
R^. There the nonlinear Schrodinger equation and complex MKdV equation play the same role as 
the MKdV equation [Ma3]. As the explicit function form of the finite type solution the nonlinear 
Schrodinger equation was obtained in [EEK], this study might be extended to that in M?. 
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Appendix 



As (2-22) was given only in improved version [BEL2] of published version [BELl], we briefly 
review its derivation here following [BEL2]. There Buchstaber, Enolskii and Ley kin started from 
a fundamental relation among the hyperelliptic sigma functions and incomplete integrals of the 
first, second and third kinds. Then they used the operator (2-28), investigated of symmetry of an 
intermediate equation and reached (2-22). 

Thus let us follow their way. First we will give the found-mental formula in the a function theory 
[Bal, BELl] 

(A-i) 



3 



where 



9 ^P, 9 .Qj 

u = ^ / du, u' = ^ / du, (A-2) 

j=l j=l 



rP,Q ^ -rP'.Q' 
-"'P',Q' — ■"'P,Q 

P /.P' /.P /.P' 



/ dui I dri + ---+ [ dug [ dvg+P^^, (A-3) 
Jq Jq' Jq Jq' 



r!i(P',Q'), 0(P,Q):=(i^-i^)|^, (A-4) 

Jq \x-xp x-xqJ 2y 



•^p 



'Q 

and Pj (Qj) is conjugate of Pj (Qj) with respect to the symmetry of hyperelliptic curve {x,y) 
{x,-y)- 

After letting all Qj to oo, we derivative it in uj, we obtain 



0(/ du + u)-Q{ du + u)+ drj~—^ n{P,q) = 0. (A-5) 



12 



Introducing R{z) = {z — Xo)F{z), for P = (a;o,yo) and Fj = {xj,yj), (A-5) becomes 



Cj{ du + u) + / dVj + Y, / 



1 ^ 

2 

i=0 



=1 " k=l ' 

Dj{F'{z))-jDj+,{F{z)) 



1 dxk 1 j/j - j/o 



i=l 



R'{z) 



Z=Xi 



(A-6) 



Then noting the fact that the left hand side is symmetrical in xq, xi, • • • Xg, while the right hand 
side is symmetrical in xi, X2 • • • Xg but does not depend on xq, (A-6) is reduced to (2-22). (2-23) 
comes from the third term in the left hand side of (A-6). 
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